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OBJECTIVE: STUDY STIFFNESS AUGMENTATION BY MATHEMATICAL DESIGN 
APPROACH: APPLY LINEAR REGULATOR THEORY WITH PROPORTIONAL FEEDBACK 


JUSTIFICATION: STIFFNESS IS READILY AVAILABLE TO DESIGNER AS PREDICTABLE 
PASSIVE CONTROL 

TIME-INVARIANT LINEAR REGULATOR— GENERAL 

SYSTEM: 

x - Ax + Bu + Dw 

CONTROLLED VARIABLES: 
y = Cx 
OBJECTIVE: 

Min J where J = 

u 

OPTIMAL CONTROL (ASSUMING w IS RANDOM;): 
u = -R -1 B T Px 
WHERE P IS SOLUTION TO 

P = -PA-A T P + PBR -1 B T P - C T QC P(t f ) = S f 

IF t f -*-<=°, GET STEADY- STATE P (AND U) FROM 
0 = -PA - A T P + PBR _1 B T P - C T QC 

POSITIVE DEFINITE P EXISTS IF 

• A IS DETECTABLE IN C, STABILIZABLE IN B 

• RESPONSE WEIGHTING MATRIX, Q, IS POSITIVE SEMIDEFINITE 

• CONTROL WEIGHTING MATRIX, R, IS POSITIVE DEFINITE 


f S f X f + 


[y T Qy + RuJ dt 



LINEAR REGULATOR ADAPTED TO STRUCTURES 


SYSTEM: 
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OBJECTIVE FUNCTION: 
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ASSUME: 


NO^a 2 ) 


• RANDOM INITIAL CONDITIONS 


• COMPLETE STATE FEEDBACK WITH NO ROTATIONAL COUPLING 

• t — *- °° (TIME-INVARIANT STRUCTURAL CHANGE) 


CONTROL: 



WHERE P 21 AND 


0 0 
-m~ 1 br _1 b t m -1 Tp 2i -m -1 br _1 b t m -1 

P 22 ARE SOLUTIONS TO 



pT A 
21 21 


- P M ^ 
21 21 21 


+ aL P T 


BR X B T (m 1 ) T P 21 + CT = 0 


( 1 ) 


AND 


P 22 A 22 + A 22 P 22 


p 22 m X BR 1 


b t (m *) 


T 


P + ( P + P 1 + 

22 v 21 21 


C 2^2 C 2^ 


= 0 


( 2 ) 


IN THESE EQUATIONS A 21 = M 1 K AND A = M 1 G 

NOTE THAT (1) IS NOT SYMMETRIC; ALSO THAT (1) IS INDEPENDENT OF (2). 
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Q WEIGHTING MATRIX CONSIDERATIONS 



• If rate and displacement considered independently and Q chosen so as not 
to couple x and x 



• For design , selection of C is governed by desired minimum response points. 
Hence, C and Q may be assigned similar functions. 

• Diagonal C and Q minimizes weighted square response at selected coordinates. 

• Choice of Q = K and Q„_ = M minimizes sum of strain and kinetic energy at 

n zz 

locations determined and (optionally) weighted by C. 
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REGULATOR FOR STRUCTURES— MODAL COORDINATES 


TRANSFORMATION x = (j)q WHERE q = q e 

K. 

WHERE 4> IS NORMALIZED = I 

AND a 1 IS ASSUMED PROPORTIONAL TO a> ± (i.E., a ± = -2^0^ OR $ T G<|> 

OBJECTIVE FUNCTION BECOMES 


J = 


NOTE THAT 



HENCE, 



WEIGHTING MATRIX Q = 




RICCATI EQUATIONS BECOME 


[r H ± M J ) 


? J n n 2 + Q 2 P 21 + P^i BR -1 -C T 1 0 2 C 1 = 0 (3) 

AND 

+ + V «- lBT? 2 2 - ( P 21 + P L + ? C 2 ) ‘ 0 (« 

WHERE P , B, R f AND C ARE MODAL EQUIVALENTS OF P, B, R, AND C. 

BY CHOOSING P, B, R, AND C DIAGONAL, WE DECOUPLE THE SOLUTION AND GET PURE "MODAL 
CONTROL." 
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CANTILEVER BEAM MODEL 



ASSUMED: 

• CONSISTENT MASS FINITE ELEMENTS 

• UNIFORM INITIAL STIFFNESS & MASS DISTRIBUTION 

• FIRST NATURAL FREQUENCY = .047 Hz (.297 rad/sec) 

PHYSICAL IMPLEMENTATION OF STIFFNESS CONTROL 
COWLING* 



COUPLING- 
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CONTROL WEIGHTING EFFECTS ON DESIGN 


UNDAMPED NATURAL FREQUENCIES 



rE : SOLUTIONS OBTAINED SEPARATELY FOR STIFFNESS AND DAMPING 

COMPARED EXACTLY TO FULL ORDER CONTROLLER SOLUTION 





STIFFNESS MATRIX COMPARISON (ASSUMED CONSTANT MASS) 
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FINAL K FOR R = .1 I 
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R WEIGHTING EFFECT ON STIFFNESS MATRIX 
(FIRST ROW ONLY SHOWN) 



RELATED SPONSORED RESEARCH 


• KAMAN AEROSPACE CORPORATION - AUTOMATED MATH MODEL 
IMPROVED FOR MATCHING EXPERIMENTAL DATA. 


• INSTITUTE FOR COMPUTER APPLICATIONS IN SCIENCE AND 
ENGINEERING - IDENTIFICATION OF EQUIVALENT PDE SYSTEMS 
TO MATCH MEASURED DATA. 


SUMMARY 


I COMPUTER PROGRAM FOR REDESIGNING STRUCTURAL MODES TO 
REDUCE RESPONSE HAS BEEN INITIATED. 

• LINEAR REGULATOR APPROACH IN MODAL COORDINATES HAS BEEN 
IMPLEMENTED. TRANSFORMATION OF SOLUTION TO PHYSICAL 
STRUCTURE IS A MAJOR PROBLEM, 

• SOLUTION OF STIFFNESS EQUATIONS AND DAMPING EQUATIONS 
CAN BE DONE SEPARATELY AS NXN SET OF (MATRIX RICCATI) 
EQUATIONS, 


PLANNED EFFORT FOR '82 


• INCLUDE MASS OF CONTROL 

• STUDY WEIGHTING TO MINIMIZE OR SELECT CROSS-TERMS 
I IMPLEMENT PHYSICAL COORDINATE SOLUTION 

I STUDY POTENTIAL FOR "BENEFICIAL" CROSS TERMS 
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